Magnetic imaging of the paramagnetic Meissner effect in the granular high-Tc 

superconductor Bi2Sr2CaCu20^ 
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We have imaged the spatial distribution of magnetic flux 
on a granular sample of the high-temperature superconductor 
Bi2Sr2CaCu2 0a: using a scanning SQUID microscope. Our re- 
sults establish the presence of spontaneous orbital magnetic 
moments which were suggested to be the origin of the para- 
magnetic response of these materials. The signature of the 
orbital magnetic moments is a rather broad distribution of 
local magnetic flelds at the surface of the sample. A simple 
model for the distribution is presented. 

Conventional superconductors generally tend to expel 
a small external magnetic field upon cooling into the su- 
perconducting state. This Meissner effect leads to com- 
plete, or (due to remnant trapped flux, e.g. in a ceramic 
sample composed of grains and voids) partial, diamag- 
netism. Therefore it came as a surprise when a paramag- 
netic piSjial was observed in ceramic Bi2Sr2CaCu202: (Bi- 
2212)ErQ. The origin of the paramamietism has been p. 
controversial subject. Braunisch et aln and Kusmartsevcl 
proposed that some form of spontaneous orbital cur- 
rents were responsible, giving rise to magnetic moments 
which could be aligned by the external field. This pro- 
posal for spontaneous orbital currents (Wohllcben Ef- 
fect) in turn led two of uaau to propose that an intrin- 
sic d^2_y2-wsLve symmetry of the superconducting state 
would naturally lead to frustrated Josephson junction 
circuits in a ceramic sarnnle where randomly oriented 
grains contact each otherEl. Although spontaneous or- 
bital currents have now been unequivocally demonstrated 
for high-temperature superconductors in controlled ge- 
ometries, and the gvid^^ice for d^^-y^-wa-ve symmetry is 
now overwhelmingtrllj, the controversy has continued. In 
part this is due to the observation of paramagnetic sig- 
nals under|-quite different conditions, e.g. in bulk Nb 
samplesll3~t3. In this letter we report the first direct 
imaging of the local magnetic flux distribution in the 
ceramics by a scanning SQUID microscope (SSM) and 
demonstrate that a polarization of the distribution of 
spontaneous fluxes is indeed responsible for the param- 
agnetic signal. 

Before proceeding to a discussion of our experiment 
we would like to remark that these two forms of para- 
magnetism in ceramic Bi-2212 and in bulk Nb sample 
can be clearly distinguished in several other ways. For 
example, the cooling rate affects the magnetic response 
differently in the two cases. Recent experimental data 
show signiflcant differences between Nb and granular 



Bi2Sr2CaCu202: (Bi2212) samples. While slow cooling 
enhances the paramagnetic signal for the granular sam- 
ple, it is diminished in the Nb sample. This clearly in- 
dicates that the equilibrium state of bath samples in a 
small magnetic field is quite differentEEl. For the Nb- 
disks Koshelev and Larkin gave an explanation based on 
the idea that during the cooling process the surface re- 
gion nucleates superconductivity before the bulk, so that 
magnetic fiux in the sample is compressed and creates an 
enhanced magnetizationEll. This compressed fiux mech- 
anism leads to a metastable state which depends on the 
cooling procedure whereas the polarization of the spon- 
taneous orbital moments is an equilibrium process. Fur- 
ther, noise measurements of the magnetization of Bi2212 
give signals which are compatible with the presence of 
spontaneous orbital momentsc3. 

We used a sample whose preparation, characteriza- 
tion and measurement of the magnetization were re- 
ported previouslyu. The naagnetic images were made 
with a high-resolution SSMc2l. This instrument uses a 
Nb-Al203-Nb low T^ SQUID fabricated on a sihcon sub- 
strate. The substrate is polished to a sharp tip spaced a 
few tens of microns from a well-shielded superconducting 
pickup loop, which is an integral part of the SQUID. The 
SQUID substrate is mounted on a flexible cantilever, ori- 
ented at a shallow angle, typically 20 degrees relative to 
the sample surface, and the sample is scanned relative to 
the SQUID, with the tip of the substrate in direct contact 
with the sample. For these measurements the sample was 
polished to a mirror flnish, and both sample and SQUID 
were coated with a thin layer to protect the SQUID sub- 
strate from abrasion. We estimate, from fits of Abrikosov 
vortices imaged using similar tip geometry and SQUID 
and sample coating techniques, that the spacing between 
the pickup loop and the surface of the superconducting 
sample is about 5/im. The SQUID signal is proportional 
to the magnetic fiux through the loop area. The present 
images were taken with a square pickup loop. In this 
geometry a single bulk Abrikosov vortex couples about 
0.5$o ("^o — h/2e) through the 8.2/im x 8.2^to area of 
the pickup loop located directly above it. 

The Fig. l(a-f) shows a series of SSM images of the 
Bi2212 sample which was cooled through the supercon- 
ducting transition temperature 84K) at different val- 
ues of an externally applied magnetic field. The images 
were taken with the field still applied and the sample and 
SQUID immersed in liquid helium at 4.2K. Each image 
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is of a square area 3mm on a side. The outlined square 
shows the 480 /Ltm x 480/zm area of the images used to 
generate the histograms of Fig. 2. We analyze only this 
area, because our SQUID has two magnetic field sensitive 
regions: the pickup loop, and the hole in the supercon.- 
ducting ground plane for the flux modulation coiSEj. 
The images were taken with the SQUID oriented verti- 
cally, with the pickup loop towards the top of the images. 
Therefore, when the pickup loop covers the region out- 
lined, the modulation hole senses areas well off the sample 
(below the bottom of the image) , and simply contributes 
a constant background signal to the image. 

In Fig.l we show the spatial distribution of the mag- 
netic flux in and around the sample. The grey contrast 
scale is chosen so that white corresponds to the largest 
and black to the smallest (often negative) flux value. In 
all cases the flux is plotted relative to the flux intro- 
duced by the external field, which sets the grey level away 
from the sample. One overall feature observable by eye 
is the difference between the paramagnetic magnetiza- 
tion (sample is brighter than the background) at weak 
fields, and a diamagnetic signal (sample is darker than 
the background), in the pickup loop at strong fields. For 
weak external fields the inhomogeneity of the magnetic 
flux is clearly visible and gives rise to a broad distribution 
of the local fluxes. 

The flux distributions relative to the external flux mea- 
sured in the outlined square are shown as histograms 
in Fig. 2. The distribution is broad, as anticipated 
above, and the average value indicates the overall re- 
sponse, which is paramagnetic for weak flelds. In Tab.l 
we list the average flux — ^ext and the standard de- 
viation — ((<i> - <&ai>)^)^^^ i^ext denotes the flux of 
the external fleld through the pickup loop). 

We flrst compare our SSM data with the-|total magne- 
tization measurement of the whole sampled. Both show 
the same qualitative dependence of the low-temperature 
magnetization on the applied fleld (see Fig. 3). In par- 
ticular, the sign of the magnetization changes for both 
measurements at the same fleld value B^xt ~ O.GG. This 
confirms that the SSM data, scanning only a part of the 
sample, is typical of the magnetization of the whole sam- 
ple. Let us now turn to the width of the flux distribution, 
i.e. the standard deviation from the average flux. The 
fleld dependence of 6^ indicates the existence of sponta- 
neous flux at low external flelds. If the flux observed were 
entirely due to flux trapped and compressed between and 
inside the grains, then we expect that both and (5$ 
would tend to zero in the zero-flcld limit. This is, how- 
ever, not the case as we illustrate in Fig. 3, where we 
observe that the zero-field limit of (5$ is finite. This can 
be readily interpreted if we assume that the fiux distri- 
bution at low fields is mostly due to spontaneous orbital 
currents for low external fields, which can flow in either 
direction. Thus we expect to see an inhomogeneous fleld 
pattern even at zero external fleld. The broadening of the 
flux distribution with increasing fleld can be understood 
as due to flux trapping and Meissner effect of the grains. 



Generally more magnetic flux concentrates in the voids 
and essentially little flux is trapped inside the grains. 
This leads to a enhancement of the contrast in the flux 
values for large external flelds and, consequently, to a 
broader distribution. 

For the low-fleld regime the flux distribution can be 
easily simulated using a model for the boundary between 
many grains. Such a grain boundary can be considered 
as a long Josephson jimction and may be described by a 
Sine-Gordon equatioiO, 

|^ = ^sin(^(a;)+0(x)), (I) 

where (p is the Josephson phase difference on the grain 
boundary and Aj is the Josephson penetration depth. 
The presence of 0- and 7r-junctions enters through 6{x), 
which assumes the values or tt as a function of the po- 
sition X along the grain boundary. This model is sim- 
ulated on a system of length L using N mesh points 
to determine ip{n) for flxed values of d{n)cB. The lo- 
cal flux <I>(n) between mesh point n and n — 1 is given by 
^(n) = $o('/5(n) — (p{n — l))/27r. The external fleld is in- 
troduced via the boundary conditions at the two ends of 
the junction {ip{N) - ip{N - l))/2n = {ip{21- ip{l))/2Tr = 
BextL/N^Q. Using a relaxation methodcJ, we calculate 
if while gradually lowering the temperature by introduc- 
ing decreasing values of Xj. In Fig. 4 we show the flux 
distribution obtained for the case of L = 100, N = 1000, 
starting with A,/ = 40, which is decreased by successive 
division by 2 to a flnal value of 0.156. The low external 
fleld is Bext = 0-l^o/Ld {d is the magnetic width of the 
junction). We get a broad distribution with a shape that 
is qualitatively similar to the experiment. The phase Lp 
has essentially a random walk like dependence on x so 
that the histogram has an approximately Gaussian dis- 
tribution. Within this simple model we can describe the 
generation of spontaneous flux and the interaction effects 
between flux lines. However, the broadening of the flux 
distribution with increasing fleld is not properly repro- 
duced because the contrast between trapped flux and the 
screening grains is not taken into account. 

In summary we would like to emphasize that the low- 
fleld data obtained by an SSM on a granular Bi2212 
sample are in very good qualitative agreement with 
the previous magnetization measurements and, in addi- 
tion, provide direct evidence for the presence of sponta- 
neous or| bi|ta l currents. On the other hand, independent 
evidenceE3'E3 indicates that the paramagnetic signal in 
Nb samples has a different mechanism, likely the one 
proposed by Koshelev and LarkinEil. It is not possible 
from this SSM measurement alone to determine the ori- 
gin of the orbital currents. However, a series of previ- 
ous experiments suggest strongly that grain boundaries 
with intrinsic 7r-phase shifts are appearing in these high- 
temperature superconductors and ar|e j^ery likely respon- 
sible for the paramagnetic responsaalla. In this sense, 
the paramagnetic response (Wohlleben effect) in granu- 
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lar Bi2212 systems is consistent in all aspects with AApf- 
planation based on dj.2_y2-waye pairing symmetryQD'EZI. 
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here. This work was supported by the Swiss National- 
fonds. 



FIG. 1. Scanning SQUID microscope images of a granular 
Bi2212 sample, cooled and imaged in various fields. Each 
image has 512x512 pixels, with 6 microns per pixel. The 
individual images are labeled by the cooling field, and by the 
maximal range of variation of the flux (in units of <l>o) in each 
case. The outlined areas are regions of the images analyzed 
further in the histograms of Fig. 2. 

FIG. 2. Histograms of the distribution of occurrence of 
each SQUID flux value in the outlined areas of Fig.l, normal- 
ized so that the integral yields unity. Each panel is labeled 
by the cooling field in Gauss. The horizontal axes correspond 
to the flux through the SQUID at a particular spot on the 
sample, relative to the background with the pickup loop far 
from the sample, i.e. "I> — $ext- 
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FIG. 3. Standard deviation of the flux distribution 
5$ = {($ - $,„)2)i/2 (circles,solid line) as a function of the 
external field. The average flux of the SSM measurement 
(diamonds,long dashed line) compares well with the appro- 
priately scaled magnetization of the whole sample (x, dashed 
line).* 



FIG. 4. Histogram for the small fleld 
Sine-Gordon model described in the text. 
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TABLE I. Average flux $av = (<!> — ^ext) measured in 
the SSM and standard deviation 5$ = {($ - ^avf)^^^ for 
different applied flelds. 



Hext 




5$/<I>o 


10 mOe 


0.0157 


0.0384 


30 mOe 


0.0128 


0.0386 


100 mOe 


0.0299 


0.0494 


300 mOe 


0.0318 


0.0675 


1 Oe 


-0.0676 


0.139 


3 Oe 


-0.51 


0.465 
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